Collective dynamics in highly relativistic heavy-ion collisions* 

A. Dumitru 

Physics Department, Yale University 
P.O. Box 208124, New Haven, CT 06520, USA 
D.H. Rischke 

bb" 

Qv^ ' RIKEN-BNL Research Center, Brookhaven National Laboratory 

a^ ■ 

Upton, NY 11973, USA 
<D ■ (February 9, 2008) 



heavy-ion collisions at CERN-SPS and BNL-RHIC energies up to = 4 GeV. 



Hydrodynamics with cylindrical symmetry in transverse direction and 
longitudinal scaling flow is employed to calculate the transverse momentum 
' spectra of various hadrons and clusters (e.g. vr, K, N, ^, A, d, He) in central 
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o 
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' We discuss the sensitivity of these spectra with respect to the initial trans- 

00 ' 

I verse density profile as well as to the choice of "freeze-out" hypersurface. 

For ^/s = 18A GeV (SPS energy) and pT < 2 GeV, overall good agreement 
^ I of the pt distributions with data is found when freeze-out occurs along the 



T = 130 MeV isotherm. Even high-py neutral pion data can be described for 



^ ' a particular choice of the initial transverse density profile. It is shown that the 



average transverse velocity {vt) of heavy hadrons and hadronic clusters is a 
good measure for the collective flow velocity. The latter is found to be rather 
similar for SPS and RHIC energies, due to the "stall" of the flow within the 
long-lived mixed phase at RHIC. In case of thermalization and hydrodynam- 
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ical expansion, the mean transverse momentum {px) increases linearly with 
the hadron mass. In contrast, the string model FRITIOF 7.02, which does 
not account for rescattering of secondary hadrons, predicts a strong depen- 
dence of {pt) on the quark composition of the hadron. Due to the different 
sensitivity to hard processes, hadrons with charm (anti-)quarks acquire signif- 
icantly more transverse momentum than hadrons without c (or even without 
s) quarks. 



I. INTRODUCTION 

A variety of experimental data for rapidity and transverse momentum spectra of various 
hadrons and clusters in relativistic heavy-ion collisions at the BNL-AGS and CERN- 
SPS indicate significant collective (longitudinal and transverse) expansion of the produced 
hot and dense fireball p-p!3|. Also at the much higher center-of-mass energy that will be 



reached at the BNL-RHIC, it is commonly expected that the parton density at midrapidity 
in central collisions is high enough to lead to rapid thermalization ||l^Hl6|. 

In the present paper we discuss the collective evolution of the system following thermal- 
ization within a hydrodynamical model. We compute the transverse momentum spectra at 
midrapidity dN/(PpTdy\y=Q, the mean transverse momenta (pr), and the mean transverse 
velocities (vt) of various hadron species for Ph + Ph at CERN-SPS energy {^/s = 18A GeV, 
entropy per net baryon s/ps ~ 40) and Au + Au at BNL-RHIC energy {^/s = 200A GeV, 
s/pb ~ 200). Finally, we also calculate the (pt) of various hadrons in Au + Au at RHIC 



energy within a string model without rescattering of produced particles (FRITIOF 7.02 ||17|| ) 
and compare its predictions to those of hydrodynamics. 
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II. MODEL DESCRIPTION 



A. Scaling Hydrodynamics 

Hydrodynamics is defined by (local) energy-momentum and net charge conservation [|18 



5^T'^^ = , d^Nt = . (1) 

j'/j-f denotes the energy-momentum tensor, and A^f the net four-current of the ith conserved 
charge. We will explicitly consider only one such conserved charge, the net baryon number. 
We implicitly assume that all other charges which are conserved on strong-interaction time 
scales, e.g. strangeness, charm, and electric charge, vanish locally. The corresponding four- 
currents are therefore identically zero, cf. eq. (|^), and the conservation equations are trivial. 
For ideal fluids, the energy-momentum tensor and the net baryon current assume the 



simple form |18 



T'^'^ = {e + p)uV-pg'^'' , N^ = Pbu'' , (2) 

where e, p, pb are energy density, pressure, and net baryon density in the local rest frame 
of the fluid, which is defined by = (p^, 0). gf^'^ = diag(+, — , — , — ) is the metric tensor, 
and = 7(1,-17) the four-velocity of the fluid {v is the three-velocity and 7 = (1 — tP)'^/"^ 
the Lorentz factor). The system of partial differential equations ([^) is closed by choosing an 
equation of state (EoS) in the form p = p{e,pB), cf. below. 

For simplicity, we assume cylindrically symmetric transverse expansion with a longitudi- 
nal scaling flow profile, Vz = z/t HTQ). At ^ = 0, equations (^ reduce to 

dtE + dT[{E + p)vT] = -{^^ + -^{E+p) , 
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dtM + dr (MvT+p) = - ( — + 7 ) M , (3) 



dtR + St {Rvt) =-[ — + -]R 



St t, 

where we defined E = T°°, M = T°^, and R = N^. In the above expressions, the index T 
refers to the transverse component of the corresponding quantity. 

The set of equations describes the evolution in the z = plane. Due to the assumption 
of longitudinal scaling, the solution at any other z can be simply obtained by a Lorentz 
boost. The above equations also imply 



dp 

dr] 



, (4) 



where t] = Artanh Vz and r = \Jt'^ — z^ . This means that on proper-time hyperbolas pressure 
gradients in rapidity direction vanish, and there is no flow between adjacent infinitesimal 
rapidity slices. However, only for net baryon free matter, ps = 0, does this automatically 
also mean that the temperature T is independent of the longitudinal fluid rapidity r]. In the 
case pb 7^ equation (H) only demands 



dT 
dr] 



+ PB 



. (5) 



dr] 

s and pb denote entropy density and baryo-chemical potential, respectively. If other charges 
like strangeness or electric charge are locally nonvanishing, additional terms appear. Equa- 
tion (j^) does not imply that the rapidity distribution of produced particles is fiat (i.e. in- 
dependent of rapidity) or that the rapidity distributions of various species of hadrons, e.g. 
pions, kaons, and nucleons, are similar. Any rapidity- dependent T and pb that satisfy eq. (|^) 
are in agreement with energy-momentum and net baryon number conservation, as well as 
with longitudinal scaling flow Vz = z/t. In this paper, however, we do not explore the ra- 



pidity dependence of the particle spectra, and thus, as discussed below, are not required to 
specify T and /x^ as functions of 77. 

The hydrodynamical equations of motion are solved on a discretized space-time grid 
{ArT = -Rt/100 = 0.06 fm. At = 0.99ArT) employing the RHLLE algorithm as described 
and tested in ||20|,^ . We have checked that the algorithm conserves total energy and baryon 
number to within 1% and that profiles of rarefaction and shock waves are reproduced accu- 
rately for various initial conditions pl|j2^ . 



B. Equation of State 

To close the system of coupled equations of hydrodynamics, an equation of state (EoS) 
has to be specified. Lattice QCD predicts a phase transition from ordinary nuclear matter to 
a so-called quark-gluon plasma (QGP) at a critical temperature of Tc = 140 — 160 MeV 
(for pb = 0). We model the high-temperature phase as an ideal gas of u, d, s quarks 
(with masses m„ = rrid = 0, = 150 MeV), and gluons, employing the well-known MIT 



bag model EoS p4|. In the low-temperature region we assume an ideal hadron gas that 
includes the well-established (strange and nonstrange) hadrons up to masses of 2 GeV (as 
listed in p5|). The phase coexistence region is constructed employing Gibbs' conditions of 
phase equilibrium. The bag parameter of B = 380 MeV/fm^ is chosen to yield the critical 
temperature Tc 160 MeV. 



C. Initial Conditions 



Equation (^) suggests that it is convenient to specify the initial conditions in scaling 
hydrodynamics on a proper-time hyperbola r = Xj. Apart from Tj, we also have to specify the 
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entropy per net baryon, s/ ps-, and the net baryon rapidity density at midrapidity, dNB/dy. 
The initial energy and net baryon densities as well as the temperature and the chemical 
potentials are then determined by the EoS. After specifying an initial transverse density 
distribution, the solution of relativistic ideal hydrodynamics is uniquely determined in the 
forward lightcone. 

For collisions at SPS energy we assume that hydrodynamic flow starts at proper time 



Tj = 1 fm/c. This is a value conventionally assumed in the literature, cf. e.g. [|T9|. Note that 
the production time of the partons in the central region might be smaller, ~ 0.2 — 0.3 fm, 
cf. e.g. However, we assume that one to two further collisions per particle are necessary 
for thermalization and onset of collective behavior. 

We further assume dNs/dy = 80, which (after resonance decays) results in ^ 60 net 
nucleons at midrapidity. If we simply scale by 1/2 (at high energy the initial isospin asym- 
metry is transferred to the pions and the nucleons become almost isospin symmetric), the 
net proton multiplicity at midrapidity agrees with preliminary NA49 data 

The specific entropy in central heavy- ion collisions at SPS energies is s/ps = 40 — 50, 
cf. e.g. refs. 0,^,^ and references therein. With this entropy per net baryon most hadron 
multiplicities as measured at CERN-SPS can be well described, e.g. the multiplicity of pions 
in central Pb + Pb at midrapidity, dN^^/dy ^ 600 

The initial energy and net baryon densities follow as = 5.3 GeV/fm^ and pb = 4.5po, 
respectively. The initial temperature and quark-chemical potentials are Tj = 216 MeV, and 
Pq = 167 MeV, ps = 0, respectively. An estimate of the initial energy density could also be 
obtained using the Bjorken formula 



dEj- 1 dy 



= ~ 3 GeV/fm^ , (6) 
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where dEx/dy = 405 GeV is the measured transverse energy in the most central Pb + Pb 
events. Since this simple formula does not account for the work performed by expansion, cf. 
e.g. it predicts an initial energy density that is almost a factor of two below that assumed 
in our calculation. We obtain dET/dy = 412 GeV on the hadronization hypersurface, and 
dEr/dy = 406 GeV on the T = 130 MeV isotherm. 

Let us now turn to the initial conditions for central Au + Au collisions at RHIC energy. 
Due to the higher parton density at midrapidity (as compared to collisions at SPS energy), 
thermalization may be reached earlier at RHIC. According to various studies |1T^- [1^] , ther- 
malization might occur within ~ 0.5 fm. Following we will assume Tj = Rt/10 = 0.6 fm 
for Au nuclei. 

Various microscopic models, e.g. PCM |1|], RQMD 1.07 [0, FRITIOF 7.02 and 
HIJING/B predict a net baryon rapidity density of dNs/dy ~ 20 — 35 and specific 
entropy of s/ ps ~ 150 — 250 in central Au + Au at y/s = 200A GeV at midrapidity. 
Perturbative QCD calculations of minijet production in leading logarithm approximation 
with a transverse momentum cutoff of po = 1 GeV yield s/pb = 150, if a soft component 
with s/ pb = 50 is assumed . 

We will employ dNB/dy = 25 and s/pb = 200. For our choice of r, and equation of 
state, the resulting energy and baryon densities at midrapidity are = 17 GeV/fm^ and 
Pi = 2.3po- The initial temperature and quark-chemical potentials follow as Tj = 300 MeV, 
/ig = 47 MeV, /is = 0, respectively. The higher specific entropy and lower net baryon 
number in Au + Au collisions at RHIC energy as compared to SPS energy enhance the pion 



multiplicity at midrapidity by ~ 50% |Q. Also, antibaryon to baryon ratios at midrapidity 
increase from ^ 0.1 to > 0.5 |33| ]. 

The total transverse energy at midrapidity on the hadronization hypersurface is 
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dET/dy = 641 GeV, and decreases to 627 GeV on the T = 130 MeV isotherm. Thus, 
the isentropic expansion reduces the energy at midrapidity by almost a factor of two 
{dET/dy{Ti) = 1153 GeV). 

In an individual event, the transverse flow velocity field can be rather complicated [|4 



however, on average over many events it vanishes at r = Xj. The radial energy and net 
baryon densities are taken to be proportional to Q{Rt — r^). Rt = 6 fm denotes the nuclear 
radius. To show the effect of initial pressure gradients on collective particle production at 
high pt we will also employ the "wounded nucleon" distribution 

e{Ti,rT) = eif (rr) , 
pBi.Ti.rT) = pif [tt) , (7) 

with 
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2\ 

ei, Pi denote the energy and net baryon densities as specified above for SPS and RHIC 
energies, respectively. 



D. Calculation of Particle Momentum Spectra 

It is generally unreasonable to assume that the fluid-dynamical description is valid in 
the whole forward lightcone. In some space-time regions the reaction rates between the 
particles are too low to maintain local thermal and/or chemical equilibrium. In hydrody- 
namical calculations one usually assumes that the transition between the fluid and the free 
streaming regime occurs within a space-time volume which is commonly approximated by a 
three-dimensional surface in four dimensional space-time ]TE|J55| . One then has to determine 
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the parameters of the free-streaming gas, characterizing its momentum-space distribution, 
from the properties of the fluid (i.e. velocity, temperature, etc.) on the assumed freeze-out 
hypersurface. From the theoretical point of view, this complicated problem is not yet satis- 
factorily solved Problems may arise, in particular, if the momentum-space distribution 



functions on the two sides of the freeze-out hypersurface (fluid and free-streaming gas) differ. 
Also, on time-like parts of the freeze-out hypersurface (i.e. with space-like normal) negative 
contributions to the spectrum of frozen-out particles may arise. For high flow velocities 
through the hypersurface, however, these contributions are negligible. 

Several attempts to understand and model the freeze-out can be found in the litera- 
ture p5|-p^. We will employ the description of Cooper and Frye |Q, which has been used 



to calculate rapidity and transverse momentum spectra at AGS and SPS energies by sev- 
eral authors p|- p!2| , p7| , p9| , ^ . For a cylindrically symmetric transverse expansion exhibiting 
longitudinal scaling flow the particle momentum spectra are given by 
aJM f 7^ f 7 1 f 7 I J ^dfT 



J d( J dcj) J drj rx T (pT cos (f)— — rriT cosh.{y — T])- 



fipX) > (9) 



d'^Prdy J J J \ dC dC 

-oo ^ 

where p^ = (m^^ cosh sin x,Pt cos Xi^^Tsinhy), and = 7T(cosh?7, frsin((/) — 
x),VTCOs{(j) — x)?sinh?7) denote the particle four-momentum and the fluid four- velocity, 
respectively. Thus, the direction of the particle momentum in the transverse plane is deter- 
mined by the angle X; while the relative angle between px and the transverse flow velocity, 
vt, is denoted by 0. C ^ [0? 1] parametrizes the hypersurface (counter-clockwise) on which 
the momentum distribution of particles is to be calculated, such that rxiC) ''"(C) specify 
the space-time points on the hypersurface. Since we focus on particles with rriT ^ Tfo, we 
approximate the phase-space distribution function / by a Boltzmann distribution, and also 



exp(-^^^^cosh(y-r7)j ^ y-^^5(y-r7)exp(^-^^^^^ . (10) 



The thermal smearing of the rapidity distribution is thus approximated by the square-root 
factor. The quahty of this approximation is less good for low-px pions which are therefore 
not considered here. 

For SPS energies, we calculate the transverse momentum spectra on the T = 130 MeV 
isotherm. This turns out to reproduce the measured spectra of vr, p, (p, A, and d reason- 
ably well, cf. below. Also, this choice is in line with freeze-out temperatures extracted by 
others BiiiHIl- 

Even without a detailed calculation of the reaction rates between various hadron species 
as a function of temperature and chemical potentials (cf. |^ for such a study and [0 for a 
first step to include it into the fluid-dynamical model) one might speculate that the freeze- 
out temperature might be higher at RHIC than at SPS. The reason is that the lower net 
baryon density has to be compensated for by an increase in the meson density, in order that 
the reaction rates be large enough to maintain local thermodynamical equilibrium. We shall 
thus calculate the particle momentum distributions at RHIC on the hypersurface where the 
fluid enters the hadronic phaseQ, A = 0. (A = ^/QGP^'ytot (jg^otes the fraction of quark-gluon 
plasma within the mixed phase.) The results are not changed if a T = 160 MeV isotherm is 
employed instead. To explore the range of uncertainty, however, we shall also show results 
for freeze-out at T = 130 MeV. 



Note that, for the EoS used here, hadronization proceeds via a shock wave [^.42]. In practice, 
we evaluate eq. (|9|) on the hypersurface corresponding to a small negative A = —0.05. This ensures 
that the fluid properties (velocity, temperature etc.) are those of hadronic matter, located either 
on or behind the shock front. 
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III. RESULTS 



A. Space-Time Evolution 
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FIG. 1. Hypersurfaces corresponding to A = 1 (boundary between pure QGP and mixed phase), 
A = (boundary between mixed phase and pure hadron phase), and the freeze-out isotherm 
T = 130 MeV; for central Pb + Pb cohisions at SPS. 

Figures |I| and ^ show the hypersurfaces where the pure QGP and the mixed phase end, 
respectively, as well as the T = 130 MeV isotherm. One observes a qualitative change in the 
space-time evolution. When comparing the evolution at SPS and RHIC, one finds that due 
to the high initial entropy density the QGP phase lives almost twice as long at RHIC than at 
SPS. Therefore, the transverse rarefaction wave can propagate further into the quark-gluon 
fiuid and accelerate it. This leads to larger transverse fiow of matter entering the mixed 
phase at RHIC than at SPS. Due to this effect the hadronization surface (A = 0) extends 
to larger r^- Note also that at SPS the fluid spends ~ 50% of its lifetime in the hadronic 
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phase, while at RHIC QGP and mixed phase occupy a much larger space-time volume. For 
comparison, we have also shown the T = 130 MeV isotherm for Au + Au reactions at RHIC. 
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FIG. 2. Hypersurfaces corresponding to A = 1 (boundary between QGP and mixed phase), and 
the freeze-out hypersurface A = (boundary between mixed phase and pure hadron phase); for 
central Au + Au collisions at RHIC. 



For our EoS, which contains all well-established hadrons up to mass 2 GeV, hadronization 
is completed at r = 1.7Rt = 10 fm (SPS) resp. r = 2.7Rt = 16 fm (RHIC). The values 



of the lifetime of the mixed phase obtained here agree with those predicted in ref. |2^ for 
the case where the ratio of degrees of freedom in the QGP and hadronic phase is 3, cf. 
Fig. 12c in [^]. If only tt or it, t], p, uj mesons are considered in the hadron gas, complete 
hadronization takes considerably longer at RHIC, r = QRt = 36 fm and r = 30 fm [Q, 
respectively (cf. also refs. P-pT],^). The back-reaction of the freeze-out on the dynamics of 
the fluid could further speed up hadronization by up to 30% [p3|. 
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B. Transverse Momentum Spectra of Hadrons and Clusters 




Pt (GeV) 



Pt (GeV) 



FIG. 3. Transverse momentum spectra of direct thermal hadrons at midrapidity (central Pb + Pb 
collisions at SPS), calculated on the T = 130 MeV isotherm. A step function has been employed 
for the initial transverse energy and baryon density profiles. 



Figure |^ compares our calculated pr spectra with h , (+) 



i, </., di,andAil 



data of NA49. The additional factor 1/2 in the A data as compared to our calculation 
approximately corrects for the S°s, which can experimentally not be separated from the As. 
The preliminary 7r° data obtained by the WA98 collaboration are also shown |^5|] (see also 
the data obtained by the CERES collaboration |^). Both the calculated spectra as well 
as the experimental data are absolute numbers without any arbitrary normalization factor. 
The reason for the underpredicted multiplicities of pions, protons, and lambdas is the fact 
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that we have not included post-freeze-out resonance decays[[ Also, we have not explored 
effects of an isospin asymmetric initial state. 

The shapes of the spectra, however, are well reproduced. In particular, one can see 
the broadening of the px distribution as the particle mass increases, and the "shoulder- 
arm" structure that develops in the spectra of the heavy clusters. This behavior emerges 



naturally in the presence of collective transverse flow |0,Q2,47,48], and has already been 



discovered experimentally at the BEVALAC and at the AGS 0. 

Note also that due to collective transverse flow the spectra at high pt are not exponentials. 
The inverse slope of the thermal pions and nucleons in the region 1 GeV< pr < 2 GeV is 
T* = 195 MeV, T;^ = 220 MeV, and increases to T* = 255 MeV, T;^ = 270 MeV, respectively, 
at higher transverse momentum, 2 GeV< < 3 GeV. This effect is even stronger in Au+Au 
at RHIC, see below. 

The apparent temperature of high-p-r pions agrees well with the preliminary result of the 
WA98 collaboration, T* = 250 MeV at pr = 2 GeV H]. This holds also for the multiplicity 
of TT^s with Pt > 2 GeV at midrapidity; while we obtain 0.32/3 = 0.11 (excluding vr^s from 
post-freeze-out resonance decays) the WA98 collaboration found 0.138 ± 0.05 (preliminary) 
TT^s in this region of momentum space. In contrast to the arguments of ref. |5^ it can thus 
not be excluded that the measured high-p^ pions can be described within a fluid-dynamical 
model. 



^The total multiplicities, including feeding from resonance decays and assuming chemical freeze- 
out at the phase boundary were calculated e.g. in refs. p J2^,3^. 
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Pt (GeV) Pt (GeV) 

FIG. 4. Transverse momentum spectra of direct thermal hadrons at midrapidity (central Pb + Pb 
collisions at SPS), calculated on the T = 130 MeV isotherm. A "wounded nucleon" distribution 
has been employed for the initial transverse energy and baryon density. 



We emphasize, however, that collective particle production at high pjn shows some sensi- 
tivity to the initial conditions. Employing a "wounded nucleon" distribution, cf. eq. (0), the 
pion spectrum (on the T = 130 MeV hypersurface) above = 1 GeV is again closer to an 
exponential, cf. Fig. ^. On the other hand, boost non- invariant initial conditions P^TT|, or 



initial transverse velocity gradients |^4|, can again enhance collective transverse expansion. 
In view of these uncertainties, we argue that collective hadron production above pt = 1 GeV 
cannot be definitely excluded. 
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FIG. 5. Transverse momentum spectra of direct thermal hadrons at midrapidity (central Au+Au 
collisions at RHIC), as calculated on the A = and T = 130 MeV hyper surf aces, respectively. A 
step function has been employed for the initial transverse energy and baryon density profiles. 



Figure ^ shows the spectra for RHIC energy as calculated on the A = and T = 
130 MeV hypersurfaces. The AA is a hypothetical bound state of two A baryons for which 
we assume a binding energy of 30 MeV. If freeze-out occurred later, e.g. on the T = 130 MeV 
isotherm like in Pb + Pb at SPS, the resulting spectra are of course "stiffer", since collective 
flow is significantly stronger than on the A = hypersurface. This leads to much broader pt 
distributions for the heavy hadrons. Observe also that the number of pions at = 4 GeV 
increases by a factor of 3. 

We would like to emphasize also another qualitative difference in the dynamics at SPS 
and RHIC, where the space-time volume of the hadronic phase might be smaller. At RHIC, 
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the outer part of the fluid (ry > 0.7Rt, cf. Fig. 0) is strongly accelerated by the rarefaction 
wave which moves into the QGP. However, once longitudinal expansion has cooled the core 
of the fluid into the phase coexistence region, the transverse rarefaction wave is stopped 
abruptly, since the velocity of sound in the mixed phase is small. Thus, if the space-time 
volume of the purely hadronic phase is small, the interior part of the fluid will remain more 
or less at rest^. We consequently expect two different slopes in the pt distributions. This is 
confirmed by Fig. |^, left panel. In the range 1 GeV< pt < 2 GeV, the apparent temperature 
of thermal pions and nucleons is T* = 220 MeV and = 233 MeV, respectively; at higher 
Pt, 2 GeV< pr < S GeV, the inverse slopes increase by roughly 100 MeV, to T* = 325 MeV, 
and = 335 MeV ! Note that the difference is much larger than in Pb + Pb at SPS, 
although the average transverse flow velocity (on the A = hypersurface) is smaller than on 
the T = 130 MeV isotherm at SPS (see below). 



^Parametrization of collective transverse flow in terms of a simple radial velocity profile [48| is 
therefore not possible on the A = hypersurface. This is in contrast to freeze-out occurring later, 
e.g. on the T = 130 MeV isotherm. 
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FIG. 6. Transverse momentum spectra of direct thermal hadrons at midrapidity (central Au+Au 
collisions at RHIC), as calculated on the A = and T = 130 MeV hyper surf aces, respectively. A 
"wounded nucleon" distribution has been employed for the initial transverse energy and baryon 
density. 



However, if the initial transverse energy and baryon density distributions exhibit smaller 
gradients, this kink in the spectra is washed out, cf. Fig. ^. For the "wounded nucleon" 
distribution, eq. (|^, the very fast tails in the velocity distribution on the hypersurface are 
absent, thus reducing the number of hadrons at > 2 GeV. Figures || and ^ show the 
possible range of collective particle production at high transverse momenta. On the other 
hand, we have found that the two initial density distributions lead to average transverse 
momenta that differ by < 10%. Therefore, we present results for (pt) only for the sharp 
density profile. 



18 



C. Mean Transverse Velocities and Momenta 
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FIG. 7. Average transverse momenta and velocities (at midrapidity) of direct thermal hadrons 
on the boundary between mixed phase and pure hadron phase, A = 0, and on the T = 130 MeV 
isotherm (for central Pb + Ph collisions at SPS). 



Figures ^ and § depict the average transverse momenta and velocities of various hadrons 
at midrapidity, as calculated from the distributions: 



.0\ dNj 
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One observes that due to the collective transverse expansion, {pt) = m{'-fTVT) increases 
linearly with the particle massQ. The {pt)tt = 470 MeV of pions at SPS is slightly over- 
predicted, NA49 measured 418 ± 30 MeV However, this is not astonishing since we 
have employed the Boltzmann approximation, and have also not included pions from post- 
freeze-out resonance decays. For all other (heavier) hadrons these effects on (pt) are less 
significant. 

On the other hand, (vt) saturates with increasing m, since the thermal velocities 



iTfo/m become negligible and {v) ~ vjiow- Thus, clusters {d, He, ...) provide the 
opportunity to determine the transverse flow velocity directly |]7|,|^, while the thermal av- 
eraging cannot be neglected for nucleons and especially for pions and kaons. To extract the 
collective flow from their (px) or (vt) requires the knowledge of the freeze-out temperature, 
and moreover might be distorted by decays of resonances l^,^ . The average transverse flow 
velocity in Pb + Pb at SPS is ~ 0.3 on the boundary between mixed and hadronic phase 
(i.e. A = 0), and increases to ~ 0.4 on the T = 130 MeV isotherm. This is consistent with 
other fluid dynamical calculations ||9|- [TT1 , |37| and numbers extracted from fits of transverse 
momentum distributions p|,^,|5|,|48[| . Our analysis does not support freeze-out temperatures 



of less than 100 MeV at SPS [pT[], at least not at midrapidity. 



^Note, however, that {'Jtvt) does in principle depend on the particle mass, since it includes flow 
as well as thermal velocity components, cf. left panel of Fig. |^. Thus, the result (pt) oc m is not 
trivial. 



20 



1.0 
0.9 
0.8 
0.7 
0.6 

A 

^0.5 
V 

0.4 
0.3 
0.2 
0.1 
0.0 





HYDRO.! =130 MeV 






HYDRO.! =130 MeV 






— ■ HYDRO. A= 






— ■ HYDRO. A= 






□ FRITIOF 7,02 






FR1!I0F 7.02 




A 

A 
\\ 








V 






■■••0— * 






/! ^ 
a" 




A 


■ 






. TT K $ 

N 


^ d AA 


He 


N AA 


He 



m (GeV) 



4 1 2 

m (GeV) 



2.5 



2.0 



1.5 > 
u 


A 

1.0 £ 



0.5 



0.0 



FIG. 8. Same as fig. but for central Au + Au collisions at RHIC. Open symbols depict the 
predictions of FRITIOF 7.02. 



Although the entropy per net baryon at RHIC is five times larger than at SPS, the trans- 
verse flow velocity increases only slightly. (In particular, the (pt) of pions almost stagnates 
from SPS to RHIC.) This is due to the fact that the transverse flow cannot increase strongly 
in the expansion of the mixed phase, since the velocity of sound is very small 
Thus, hydrodynamics predicts that the transverse flow velocity at RHIC is comparable to 
that at SPS, if the EoS exhibits a first-order phase transition. If freeze-out at RHIC occurs 
as early as on the A = hypersurface, the transverse flow velocity could even be smaller 
than at SPS, cf. Figs. and § 

It is also interesting to note that both at SPS and RHIC a significant part of the transverse 
flow is already produced within the QGP and (to a smaller part also within) the mixed phase, 
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compare the dashed curves in figs. ^ and || with the sohd fines. Tfiis befiavior is different 
to tfiat observed witfiin tfie RQMD 2.3 model, wfiere tfie transverse fiow of baryons emerges 
mainly from tfie late stage of tfie reaction |TH]. Tfie reason is tfiat witfiin RQMD 2.3 tfie 



fiadron gas is not in kinetic equilibrium at times t < 5 fm, wfiile we assume early kinetic 
equilibrium of tfie QGP in tfie central region, Tj = 1 fm. 

Tfie flow of tfie (f) meson is closely related to tfiis point. Even if mesons do not rescatter 
witfi surrounding fiadronic or mixed pfiase matter, tfieir {pr)-, (vt) or inverse slope in AA 
collisions sfiould be affected by collective flow if a QGP witfi s quarks in local tfiermal 
equilibrium is created. Tfiis is due to tfie fact tfiat tfie transverse momenta of s quarks are 
subject to a Doppler sfiift according to tfie collective flow velocity built up in tfie QGP pfiase. 
According to Fig. |^, (pt)^ ~ 680 MeV and (fT)^ ~ 0.5 at SPS, even if tfie 0's decouple from 
collective flow immediately at tfie boundary between mixed and fiadronic pfiase. Indeed, 
preliminary NA49 data sfiow tfiat tfie rougfily fits into tfie pr(^) systematics predicted 
by fiydrodynamics, as already indicated by Fig. |^. Note tfiat tfie "conventional" mecfianism 
to generate pt in cascade models, i.e. resonance decay or color flux-tube fragmentation 
followed by additional rescattering witfi tfie measured free (or additive quark model) cross 
section, apparently works for tt, K, N, A but not for tfie (f) (cf. also tfie discussion in 
ref. 0). 

For comparison, we fiave also calculated tfie average transverse momentum and ve- 
locity of direct fiadronsQ (i.e. resonance decays were also switcfied off) at RHIC witfiin 
FRITIOF 7.02 cf. flg. ^ In tfiis string model, secondary fiadrons do not rescatter. Tfius, 
collective transverse flow does not develop. One observes tfiat even tfiougfi fiard scatterings 



^Clusters {d, AA, and He) are not computed within this model. 
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are taken into account, the {pr) of N , (p, A, and Q is not higher than in case of hydro dynamic 
expansion with freeze-out on the A = hypersurface. This difference would be even bigger 
if collective flow set in earlier (remember that we assume Xj = 0.6 fm). 

Only the (pr) of charmed hadrons, D meson and Ac (2285) baryon^, is strongly enhanced 
by (semi-)hard scatterings, and is close to the solid line in fig. |^. These very heavy hadrons 
are produced by fragmentation of semi-hard jets and acquire high transverse momentum. 
Thus, in FRITIOF 7.02 the (prp) of various hadrons does not increase linearly with mass but 
depends strongly on the sensitivity to (semi-)hard scatterings, and thus to the flavor content. 
In contrast, if D mesons thermalize [0] and follow the collective transverse flow, their (pt) 
is "pulled" back onto the linear {pt){''^) systematics. 

IV. SUMMARY 

In this paper we have investigated the collective evolution of the central rapidity region 
with initial conditions as expected for heavy-ion collisons at CERN-SPS and BNL-RHIC 
energies. Even in the case where the phase transition (assumed to be of first order) from the 
plasma of m, d, s quarks and gluons to the hadron gas proceeds in local thermodynamical 
equilibrium, collective transverse expansion and equilibration of heavy mesons and (anti- 
)baryons reduces the hadronization time to < 3-Rt (for Au nuclei). Nevertheless, at RHIC, 
the lifetimes of the QGP and the mixed phase are significantly longer than at SPS. However, 
at SPS the midrapidity region evolves further into the hadronic phase, which occupies a 



Pions and kaons are not very sensitive to the dynamical evolution and are left out of the discussion 
for the moment. 



23 



larger space-time volume than the QGP or mixed phase. As a consequence, the transverse 
flow velocity and the mean transverse momentum of the particles increase substantially 
during this stage of the reaction. Without this "hadronic afterburner" the hydrodynamical 
model as discussed in this paper would not describe the measured single-particle transverse 
momentum distributions. 

At RHIC, one might expect a higher freeze-out temperature than at SPS, since the baryon 
density is smaller. If, for example, freeze-out occurs directly after hadronization is complete 
(i.e. on the A = hyper surf ace), the average transverse flow velocity turns out to be similar 
or even smaller than at SPS. 

We have calculated the transverse momentum distributions of various hadrons up to 
Pt = 4: GeV on the T = 130 MeV isotherm and on the A = hypersurface for SPS and RHIC, 
respectively. Large gradients in the initial transverse density distributions (e.g. a sharp step 
function) can strongly enhance particle production at high pt- The recently observed tt^ 
distribution (up to pt = 4 GeV) in Pb + Pb at a/s = ISA GeV can in principle be described 
by a hydrodynamic expansion. However, the high-pr tails of the hadron distributions (above 
Pt — 2 GeV) are sensitive to the initial pressure gradients. 

Within hydrodynamics the mean transverse momentum of hadrons increases linearly 
with the particle mass and is independent of other quantum numbers (provided the cor- 
responding hadron species is in local thermal equilibrium). Heavy hadrons and clusters, 
if identified, allow for a direct measurement of the collective transverse expansion velocity, 
without requesting knowledge of freeze-out temperatures. 

Finally, we have compared the fluid-dynamical results to those obtained within the string 
model FRITIOF 7.02. Due to the increasing "sensibility" on (semi-)hard processes, the (pt) of 
hadrons with c quarks increases more strongly than the {pr) of hadrons with only light (anti- 
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)quarks. Consequently, there is a strongly nonlinear relation between the mean transverse 
momentum and the particle mass within this model. 
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